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Abstract: We consider two entangled accelerating qubits coupled with real scalar helds, 
each described by the Unruh-Wald model. It is demonstrated that because of the Unruh 
effect of the helds, the bipartite entanglement between the two qubits suddenly dies when 
the acceleration of one or more qubits are large enough. We also consider three entan¬ 
gled accelerating qubits in GHZ state and in W state, with equal acceleration-frequency 
ratio, and found that in either state, the tripartite entanglement suddenly dies at a certain 
value of acceleration-frequency ratio. The equivalence between the Rindler metric and the 
Schwarzchild metric in the vicinity of the horizon of a black hole implies that for two en¬ 
tangled qubits outside a black hole, the entanglement suddenly dies when one or both of 
the qubits are close enough to the horizon, while the three entangled qubits in GHZ or W 
state, the tripartite entanglement suddenly dies when these qubits are close enough to the 
horizon. 
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1 Introduction 

Recent two decades witnessed both intensive and extensive investigations on quantum en¬ 
tanglement, “the characteristic trait of quantum mechanics” in the words of Schrodinger [1]. 
More recently, this trend has been extended to the realms of high energy physics. The rela¬ 
tivistic effects of acceleration and gravitation on quantum entanglement were investigated, 
shedding new light on the subject of quantum effects of gravity. 

A particularly concerned subject is the Unruh effect, that is, the particle content of a 
field is observer-dependent, and thus an accelerating detector in the Minkowski vacuum of 
a field feels a thermal bath of particles of this field [2-8]. The consequences of the Unruh 
effect on various kinds of entanglement have been studied. For example, because of the 
Unruh effect, a state of two field modes that is maximally entangled in an inertial frame 
becomes less entangled when one of the modes is observed by an accelerating detector, 
and degrades with the increase of acceleration, towards zero at the limit of infinite accel¬ 
eration [9]. Analogous entanglement degradation occurs for field modes observed by two 
detectors close to the horizon of a black hole, with one of them freely falling while the 
other barely escapes [9, 10]. This problem was then extended to the case of three entan¬ 
gled modes, and it was found that when one of the modes is observed by an accelerating 
detector, the tripartite entanglement, which cannot be reduced to all kinds of bipartite en¬ 
tanglement, does not approach zero in the infinite acceleration limit [11, 12]. Extension of 
such studies to Fermion fields was also made [13-16]. These results are consistent with the 
general feature that entanglement of a field state depends on the choice of single particle 
modes [17]. Entanglement between detectors have also been studied in the case that one 
of the two entangled detectors accelerates while the other moves uniformly, and was found 
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to exhibit entanglement sudden death [18]. One approach is to model the detectors as har¬ 
monic oscillators [19], for which calculations were also made on the physical details [20, 21]. 
Another approach is to model the detectors as qubits [22]. Yet another approach is to 
consider the detectors as an open quantum system [23-25]. Other quantum informational 
quantities such as discord [26, 27] and quantum Fisher information [28, 29] have also been 
studied. So far there was no work on the case that more than one of the entangled detectors 
accelerate. 

In this paper, we consider two causally separated but quantum-entangled qubits, each 
of which independently accelerates and is coupled with a scalar held as described by the 
Unruh-Wald model. Our work originated in generalizing a previous work on the decoherence 
of one qubit due to acceleration [30]. We show that because of the Unruh effect of the helds, 
the entanglement between the qubits vanishes at hnite values of acceleration instead of in 
the inhnite limit, i.e. exhibit entanglement sudden death. Our result implies that outside 
of the black hole, the entanglement between the two entangled detecting qubits vanishes 
when one or both of the qubits are close enough to the horizon. We also report a result on 
the tripartite entanglement in three entangled accelerating qubits. 

2 Formalism 

Let us consider two qubits A and B far away from each other, that is, we assume there is 
no causal connection between the two qubits. For each qubit, we apply the model of the 
Unruh and Wald [6]. 

The Hamiltionian of each qubit q (q=A,B) is 

Hg = ngQlQg, (2.1) 

where the creation operator Qq and annihilation operator Qq are dehned by Qq\0)q = 
Q\\l)q = 0, Qq|0)g = |l)q and Qq\l)q = |0)q, with Subscript q=a,b. Ug gives the energy 
difference between eigenstates |l)g and |0)g. 

Each qubit is locally coupled with a held $g within a small region around it, the 
interaction Hamiltonian being 

Hlg{tq) = €q{tq) [ $g(x) [V’g (x)Qg -|-V’g (x)QJ] (2.2) 

7 Eg 

where x and tq are spacetime coordinates in the comoving frame of the qubit, the integral 
is over the spacelike Cauchy surface Eg at given time tq, eq{tq) is the coupling constant with 
a hnite duration of qubit-held interaction, 'ipq{'x) is a smooth function nonvanishing within 
a small volume around the qubit [5]. The helds and could be the same or different. 

We presume that the distance between the two qubits is so large that there is no 
physical coupling or inhuence between the neighboring helds of the two qubits during the 
interaction times. The total Hamiltonian is simply Ha + + Hb + + Hjg, 

where is the Klein-Gordon Hamiltonian for $g. In the Minkowski spacetime, each qubit 
is conhned in its own Rindler wedge and possesses boost Killing helds which are timelike. 
The only extra constraint on the trajectories of the two qubits is that the time interval of 
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the interaction between each qubit and its neighboring field multiplied by the speed of light 
is smaller than the shortest distance between the interaction regions of the two qubits. 

Therefore after a time duration longer than the interacting times Tq S> the state 

of the whole system in the interaction picture is transformed by 

Ua <8) Ub, 

where Uq is the unitary transformation acting on qubit q and the field ^q in its neighboring 
region, as given by the Unruh-Wald model. To the first order [6], 

Uq^l-i j d>go(i',x)eq(t')[<3goe"**^''*Vg(x) + 'tp*q{yi)]^/^d?xdt', (2.3) 

where d>qo and QqQ are ‘hg and Qq for Cq = 0, i.e. when the qubit-field coupling is turned 
off. It can be obtained that 

Uq^l + iQqoa\r*q) - iQloa{r*q), (2.4) 

where a(rp and a^(rp are the annihilation and the creation operators of T*, with 

Tq{x) = -2i j[G r{x-,x') - Ga{x-, x')]eq{t')e''^^'tjjq{x.')y/^d‘^x', (2.5) 

Gnq and Gaq being the retarded and advanced Green functions of the field <I>g, respectively. 
For any mode we can write 

a(r*) = {T*q,x,)a{x,) + {T*q,x'M^,). (2.6) 

at(rp = (q,xg)*at(xg) + (r*,x;)*at(4), (2.7) 

where x'q is some mode orthogonal to Xq- | q)Xq]dS^^, where 

Eg is some Cauchy surface. This inner product can be assumed to be negligible unless Xg 
is at a frequency ~ fig. Therefore, each qubit q is only coupled with the held mode xi^q) 
with frequency fig, which is further assumed to be nondegenerate. Hence we only need to 
consider xOa X^Ib i^ studying the qubits A and B. All the other modes are decoupled 
with the qubits. 

Now we consider the Fock state |n)o^, containing n particles in the mode xi^q) of th® 
held $g, as observed in the Rindler wedge conhning qubit q, 

aiT*q)\n)n^ = fiqa{nq)\n)n^ = Hqy /n\n - l)n^, 

= //gaf(Hg)|n)o^ = //*Vra + l|n l)o,, 

where /ig = (r*,xo,) = /eg(f)e*^-jVg(x)xo,(Tx)v^d'^a;. 

Hence Uq evolves only the qubit q and the mode Xo,; while the other modes of ‘hg are 
not affected and can be ignored, 

Un,\0)q\n)n^ = |0)g|n)Q^ - i^/^Hq\l)q\n - l)n,, (2.9) 

Un,\l)g\n)n, = \i-)q\n)n, + Wn + l/ig|0)g|n + 1 )q,. (2.10) 
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3 Entangled states of the two detecting qubits 

We now suppose the initial state of the two qubits to be 


l^i) = a|0)A|l)B + /3|l)A|0)B, (3.1) 

where a and /3 are superposition coefficients satisfying |ap + |/3p = 1. The results for the 
initial state of the form of a|0)yi|0)B + P\^)a\^)b are similar. Without causal connection 
between the two qubits or between the fields, each qubit detects a thermal bath of the 
Unruh particles determined by its own acceleration. With each qubit in its own Rindler 
wedge, the initial state of the whole system, as observed by the observers comoving with 
the qubits, is described by the density matrix 

Pi = \^i){^i\® P^A® P\ (3-2) 

where 

pn, = C, (3.3) 

riq 

ttq is the acceleration of qubit q, Cq = \/l — p' jg ^-j^g state of the other decoupled 

modes and is ignored henceforth. 

The final state of the system in the interaction picture is 

Pf = uIu\p^UaUb, (3.4) 

which can be evaluated by substituting Eqs.(2.9) and (2.10). Subsequently by tracing out 
the fields, we obtain the reduced density matrix of the two qubits, with respect to the 
comoving observers 


PAB — C'IC'b 

nA,nB 


ZnA,nB 


x{[\a\‘^{nB + 1)\pb\‘^ + \/3\‘^{nA + l)|/iAp]|00)(00| 
+a/3*|01)(10| + /3a*|10)(01| 


+ [\a\‘^ + \l3\^nB{nA + l)|AtBp|liAp]|01)(01| 
+ [\a\‘^nA{nB + 1)\pa\^\pb? + |/3p]|10)(10| 
+ [\a\‘^nA\pA? + |/3pnB|//Bp]|ll)(ll|}, 


where = 1 + nA\pA? + nB\pB? + nAnB\pA?\pB? + \a\‘^{\pB? + nA\pA?\pB?) + 

\P?{\PA? + nB\pB?\pA?). 

In the case that qubit A moves uniformly while qubit B accelerates, pab is 

„2 

PAB - ®b1^ 

riB ® 

x{[\aWnB + l)\pB\‘^ + W\pA?]mm 

+a,d*|01)(10| + /3a*|10)(01| 

+ [|ap + |/ 3 preB|/xsp|//Ap]| 01 )( 01 | 

+|/ 3 pll 0 )( 10 | + |/ 3 pnB|/XBnil)(ll|}, 


where = 1 + nB\pB? + \a\‘^\pB? + \l5\^\pA? + \l3\^nB\pA?\pB?■ 

We now study the correlation and entanglement in pab- Note that the entanglement 
and correlation are respectively the same in Schrodinger and interaction pictures. 
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4 von Neumann entropy S{pab) and mutual information I{A\ B) 

The von Neumann entropy 

S{pab) ^ -ArpAB'^OgpAB (4.1) 

is a measure of mixture of pAB- On the other hand, in the Minkowski frame, the state of 
the whole system is a pure state, and S{pab) quantifies the entanglement between the two 
qubits on one hand, and the fields on the other. 

From S{pab), we also calculate the mutual information 

I{A:B) = S{pa) + S{pb)-S{pab), (4.2) 

where 

PA = TrBpAS, (4.3) 

PB = TtaPAB- (4.4) 

I {A : B) is the difference between the sum of the entropies of A and B as a whole on 
one hand, and the entropy of A plus B as a whole on the other, and is thus a quantification 
of the total correlation contributed by both entanglement and classical correlation. In the 
numerical calculations throughout this paper, the bases of the logarithms are chosen to be 
2, and the parameters pA, Pb are both set to be 0.1. 

S{pab) and I{A : B) are depicted together in Figs. 1 to 3 for three cases. The result 
for the case of qubit A uniformly moving while qubit B accelerating is shown in Fig. 1. 
The special case of two maximally entangled qubits with one of them uniformly moving 
was previously studied by using a different approach and assuming no coupling between the 
uniformly moving qubit and the held [7]. In our studies, both qubits always couple with the 
helds. Fig. 2 depicts the result for the case that the acceleration-frequency ratios of the two 
qubits are always equal. Fig. 3 gives the results for various given values of qa/^a- As shown 
in these hgures, when both ua/^a and ub/^b are close to 0, pab is close to the pure state 
iTj), hence pA is close to |ap|0)(0| + |/3p|l)(l| while ps is close to |/3p|0)(0| -|- |ap|l)(l| , 
therefore S{pab) is close to 0 while I{A : B) is close to —2 |q:P log |ap — 2|/3p log [/Ip. With 
the increase of one or both of the acceleration-frequency ratios, S{pab) quickly increases up 
to a maximum while I {A : B) quickly decreases down to a minimum. The actual extreme 
values depend on the details of the dynamics, but we can make the following estimation. 
If PAB is maximally mixed, S{pab) reaches the absolute maximum 2 while S{pa) and 
S{pb) reach the absolute maximum 1, because pab is 4-dimensional while pA and ps are 
2-dimensional. Therefore the minimal value oi I {A : B) is near 0. When the acceleration- 
frequencies further increase, S{pab) slowly decreases while I {A : B) slowly increases. In 
the limit of as/^B —t oo while ua/^a = 0, pab —t |ap|00)(00| -|- |/3p|ll)(ll|. In the limit 
of both ua/^a and aB/^B approach oo, pab —t |ap|I0)(I0| -|- |/3p|0I)(0I|. In both of these 
two limits, S{pab), S{pa), S{pb) and I{A : B) all approach —|aplog|ap — |/3plog|/3p. 

5 Entanglement between qubits A and B 

Now we turn to the entanglement between qubits A and B, which are in the mixed state 
Pab- For the two-qubit mixed state pab, a measure of the entanglement is the logarithmic 
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Figure 1 . S{pab) [plots (1) to (4)] and mutual information I{A ■. B) [plots (5) to ( 8 )[ as functions 
of the acceleration-frequency ratio of qubit B, in the case that qubit A moves uniformly, for different 
initial states. (1,5): a = (2,6): a = 0.4, (3,7): a = 0.2, (4,8): a = 0.1. 



Figure 2 . S{pab) [plots (1) to (4)[ and mutual information I{A : B) [plots (5) to ( 8 )[ as functions 
of the acceleration-frequency ratio, which is assumed to be the same for the two qubits, for different 
initial states. (1,5): a = l/-\/2, (2,6): a = 0.4, (3,7): a = 0.2, (4,8): a = 0.1. The extreme values 
are insensitive to the initial state. 

negativity 

^ogWpAsW’ ( 5 - 1 ) 

where ||/ 0 ^el| is the trace norm of p^, which is the partial transpose of pab [31, 32]. 
Another entanglement measure is the concurrence 

C{p) = max{0, Ai - A 2 - A 3 - A 4 }, (5.2) 

where Aj (i = 1, 2 ,3,4) are decreasingly ordered eigenvalues of the matrix 

with pab = (o'y (S> ^y)PAB(^y ® being the y-component Pauli matrix [33]. 
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Figure 3. S{pab) [plots (1) to (4)] and mutual information I{A : B) [plots (5) to (8)[ as functions 
of the acceleration-frequency ratio of qubit B, for different given values of acceleration-frequency 
ratio of qubit A: (1,5): ua/^^a = 0, {2,Q):aA/i^A = 100, (3,7): ga/^a = 200, (4,8): ua/^a = 300. 
It is set that a = [3 = 



Figure 4. Logarithmic negativity [plots (1) to (4)[ and concurrence [plots (5) to (8)[ as functions of 
the acceleration-frequency ratio of qubit B, in the case that qubit A moves uniformly, for different 
initial state. (1,5): a = l/-\/2, (2,6): a = 0.4, (3,7): a = 0.2, (4,8): a = 0.1. 


As shown in Figs. 4 to 7, the logarithmic negativity and concurrence decrease with 
the increase of the acceleration-frequency ratio of each qubit and, especially, suddenly dies 
at a finite value of the acceleration-frequency ratio. As can be seen in these figures, the 
acceleration-frequency ratio of one qubit at which the entanglement suddenly dies decreases 
with the increase of that of the other qubit. When one of them is zero, the other must be 
larger than some finite value. In Figs. 6 and (7), the 3D plots of logarithmic negativity and 
concurrence are symmetric with respective to the plane = ub/^b, as can be seen 

in the the expression of pAB- These 3D plots also indicate that entanglement sudden death 
occurs on a curve of aA/^A and ub/^b- 
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Figure 5. Logarithmic negativity [plots (1) to (4)] and concurrence [plots (5) to (8)[ as functions 
of the acceleration-frequency ratios of A and B, which are assumed to be the same, for different 
initial state. (1,5): a = (2,6): a = 0.4, (3,7): a = 0.2, (4,8): a = 0.1. 



Figure 6. Logarithmic negativity as a function of qa/^a and qb/^b- a = (3 = 


6 Tripartite entanglement in three entangled accelerating qubits 

We have also extended our study to three accelerating qubits A, B and C, by using the 
formalism similar to the two-qubit case above. It is well known that there are two types of 
3-qubit states, GHZ state 


1 


and W state 


\GH Z) ABC = ^{\000) + \lll)), 


\W)abc = ^(1001) + 1010) + |100)), 


( 6 . 1 ) 

( 6 . 2 ) 


each representing a different type of tripartite entanglement [34]. 

We apply the above formalism to three entangled qubits A, B and C, as described by 
Eq. (2.1), each of which is locally coupled with a field (q=A, B, C). For the reason 







Figure 7. Concurrence as a function of the acceleration-frequency ratios of A and B. a 


/3 


1 

C 2 - 


given above, only the mode Xfig needs to be considered. In this way, for GHZ state, we 
obtain the density matrix of the qubits. 


PABC {GHZ) 


1 

2 


/^2 /o 2 /^2 


E 

nA,nB,nc 


^-2-K{nA^Al a-A+riB^B I a-B+ric^C I a-c) 

ZnAriBnc 


[(1 + {nA + 1) {riB + 1) (nc7 + 1) \piA?\pB?\pc?) |000) (000| 

+ (1 + nAnBnc\pA?\pB?\pc?) |111) (111| + |111) (000| 

+ | 000 ) ( 111 | + {nA\pA? + (ns + 1 ) {nc + 1 ) \pb?\pc?) | 100 ) ( 100 | 
+ (nsl/r-Bp + {ba + 1) (nc + 1) |/^.4p|/^cp) |010) (010| 

+ (ncl/^cP + (nA + 1) (ns + 1) |/^a|^|/^bP) |001) (001| 

-h (nAnsIlXAl^l/isP + (nc + 1) \pc?) |110) (110| 

+ (nAnc|/^Ap|/^cP + (ns + 1) l/^Bp) |101) (101| 

-h (nsncl/isl^llicl^ + (nA + 1) \PA?‘) |011) (Olll] , 


(6.3) 


where 


^riAriBnc —2 + (2nA + l)|/J'Ap + (2ns + ^)\pb\'^ + (2nc + l)|/^cP 

-h (2nAns + nA + ns + 1)\pa\‘^\pb\‘^ + {‘^nAUc + nA + nc + l)|AiA|^|lic|^ 

-I- ( 2 nsnc + ns + nc + l)|/is|^|lic|^ 

-h {InAUBnc + riABB + nsnc + riAnc + nA + ns + nc + 1)\pa?‘\pb\^\pc\^■ 

(6.4) 


If the three qubits are in W state, their density matrix can be obtained as 
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2-K{nA^A/CLA+riB^B / CLB+nc^C / a-c) 


PABC{W) 


_/^2 /^2 /^2 \ '' 

UA^riB^nc 


z. 


[| 001 )( 010 | 


UAriBnc 


+ |001)(100| + |100)(001| + |100)(010| + |010)(001| + |010)(100| 

+ nB\pB\^\On){lW\+nA\pA\Aoi){nO\+nc\pc\‘^\OU){Wl\ 

+ nA\pA\Aw){Wl\+nc\pc\Ao^){On\+nB\pB\Aw){On\ 

+ {{'iT'A + 1) l/^Ap + (bb + 1) + (nc + 1) l/^d^) |000)(000| 

+ (l + {ha + l)nc|//ApdcP + {bb + 1 ) bc \ pb \‘^\ pc \‘^) |001)(001| 

+ (l + {ha + 1) bb\pa\‘^\pb\‘^ + bb (bc + 1) |010)(010| 

+ {{uA + 1 ) bbbc\pa\‘^\pb\‘^\pc\‘^ + bbIpb]'^ + rec |/ icP ) | 011 )( 011 | 

+ (l + UA {bb + 1 ) + ba {bc + 1 ) dApdcP) | 100 )( 100 | 

+ {ba {bb + 1) bcIpaI'^Ipb^Ipc]^ + ba\pa\'^ + Bc\pcf) |101)(101| 

+ {baIpa^ + bbIpb^ + BABB {bc + 1) IpaI'^Ipb^IpcA |110)(110| 

+ {babb\pa\‘^\pb\'^ + babcIpaI'^IpcI'^ + bbbqIpbI'^IpcI'^) | 111 )( 111 |] , 

(6.5) 


where 

ZnAnsnc =3 + {^Ba + l)dAp + {^Bb + l)|h's|^ + (371(7 + l)|/i(7p 

+ {Suabb + ba + bb)\pa\‘^\pb\‘^ + {Sbbbc + bb + nc)\pB\‘^\pc\‘^ 
+ {Shabc + ba + bc)\pa\‘^\pc\‘^ 

+ {Suabbbc + BABB + bbbc + babc)\pa\‘^\pb\‘^\pc\‘^ ■ 


The tripartite entanglement is the genuine three-party entanglement that cannot be 
reduced to any bipartite entanglement. We use the negativity three-tangle as the measure 
of the tripartite entanglement [35], which is dehned as 



TT = + 7r_B -1- TTc), 

(6.7) 

where 

'^A = J^A{BC) ~ ^AB ~ J^ACi 

(6.8) 

with 

^A{BC) = Wp^ABcW ~ 

(6.9) 


Ads = Wp^abW ~ 

(6.10) 


Here pABC is the density matrix of the three qubits, pAB is the reduced density matrix of 
A and B, other quantities are similarly defined. 

For simplicity, here we only present the result for the case that the acceleration- 
frequency ratios of the three qubits are the same. As shown in Fig. 8. For either GHZ or W 
state, the negativity three-tangle decreases with the increase of the acceleration-frequency 
ratio, and suddenly dies at a certain value. We have also found that entanglement sud¬ 
den death generally occurs when at least two qubits accelerate, no matter whether the 
accelerations are equal. More details will be discussed elsewhere. 
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100 200 300 400 

a./n. = ajn„ = a 


Figure 8. Negativity three-tangle as a function of the acceleration-frequency ratio of the three 
qubits in GHZ state (1) and in W state (2). 


7 Discussion and summary 


Now we come to the black hole. The spacetime outside its horizon is described by the 
Schwarzschild metric 

= (1 - —)dt‘^ - (1 - - ydO"^ - r^sin^ddi/)^, (7.1) 

r r 

where the notations are standard. The proper acceleration of a static observer at r is 


a = 




(7.2) 


It is well known that near the horizon r = 2m, the Schwarzschild metric can be approxi¬ 
mated as the Rindler metric [10, 36]. The closer r is to the horizon, the larger the acceler¬ 
ation. The uniform movement corresponds to free falling into the black hole. 

Consider the entangled states studied above. Suppose at most one of them freely fall 
into the black hole, while the other is near the horizon r = 2m. According to the calculation 
above, we know that the entanglement between the qubits suddenly dies when one or more 
accelerating qubits are close enough to the horizon. For each qubit near the horizon, we 
have 


r ~ 2 m[l 


1 

(4ma)2 


-1 


(7.3) 


from which the location of the entanglement sudden death can be determined. 

Finally, one may wonder the reason of the entanglement sudden death as studied here. 
We think the Bosonic fields act as a drain of the entanglement originally exists between the 
qubits, because there are infinite number of Fock states |n) for each Bosonic field mode. 
In contrast, we conjecture that there is no entanglement sudden death if the fields are 
Fermionic as there are only two Fock states |0) and |1) for each Fermionic field mode. The 
absence of entanglement sudden death was recently noted in entanglement between Unruh 
modes [16]. Our work implicates that entanglement sudden death of the qubits can act 
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as a probe of the nature of ambient fields coupled with the qubits and the existence of 
acceleration or gravity. 

To summarize, we have studied the entanglement of two accelerating qubits coupled 
with scalar fields, and demonstrate the occurrence of its sudden death. We also found the 
entanglement sudden death of tripartite entanglement of three accelerating qubits in GHZ 
and W states. These results imply the entanglement sudden death of field-coupled qubits 
near the horizon of a black hole. This work might be useful to the issue of black hole 
hrewall [37] or energetic curtain [38]. 
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